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The phase-integral method (PIM) is an asymptotic method of the geometrical optics or semi-classical type for solving 
approximately, but in many cases very accurately, a wide class of differential equations in physics. Unlike the related 
(J)WKB method, the higher-order corrections in the PIM can be generated from a generic, unspecified base function, 
providing added symmetry and flexibility. However, with the conventional approach of using the next-to-lowest (third) 
order correction to the integrand in the phase integral as a platform for calculating higher (fifth, seventh, ninth,. . . ) order 
corrections, the higher-order calculations very often become quite complicated. 

We therefore introduce a new platform function, which considerably simplifies the calculation of the third-order 
contribution for a wide range of problems. We also present directly integrable conditions for the phase integral, which 
so far seem to have gone unnoticed. 

For a large number of observables, our analysis makes a clearer distinction between physical and, in a sense, unphys- 
ical contributions. 

PACS numbers: 03.65. Sq, 02.30.Hq, 02.30.Mv, 02.60.Nm 
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I. INTRODUCTION 



In their pioneering 1936 article on the calculation of 
Coulomb wave functions in non-relativistic quantum mechan- 
ics, Yost, Wheeler and Breit 1 commented on the fact that the 
WKB method in the first-order approximation sometimes has 
an unsatisfactory accuracy. Moreover, they found that the re- 
placement, originally suggested by Kramers, 2 of /(/+ 1) by 
(/ + 1 /2) 2 where / is the orbital angular momentum quantum 
number, does not always improve the results obtained in the 
first-order WKB approximation. This situation can mostly 
be remedied by going to the next-order approximation. The 
phase-integral method of Froman and Froman 3 provides a 
straightforward structure for calculating the contributions of 
order five and higher. For this, Ref 4 introduces a platform 
function £o(z) which, apart from a factor 1/2, is the third- 
order contribution. 



In the Froman and Froman PIM, the third-order contribu- 
tion is the most important correction. Hence, it is desirable to 
find ways to simplify the calculation of this contribution. For 
this purpose, we here introduce a more flexible platform func- 
tion P.s(z) which generalizes £o(z) in a way that depends on 
the number of singularities in the integrand and which has the 
potential to simplify the evaluation of higher-order correction 
terms. As a result, the choice of the so called base function, 
usually denoted Q(z), becomes more transparent. 



II. HIGHER-ORDER PHASE-INTEGRAL 
APPROXIMATIONS AND THE PLATFORM FUNCTION 



Let us consider the differential equation 



-^r+R(z)y = o 



(i) 



where R (z) is a meromorphic function of the complex variable 
Z. Following Ref. 4, we introduce into equation (1) a "small" 
book-keping parameter X, ultimately to be set equal to unity, 
such that we obtain the auxiliary differential equation 



dV 
di 2 " 



Q 2 {z) 
X 2 



R(z)-Q z (z))w = o 



(2) 



which goes over into equation (1) when X — I. The meromor- 
phic function Q(z) is the unspecified base which determines 
the phase-integral approximation of order one. The auxiliary 
differential equation (2) has two linearly independent solu- 
tions 



U = q~Hz) 
f_=q-i( z )e 



,+iw(z) 
-iw(z) 



where 



w(z) = J q(z)dz 



(3a) 
(3b) 



(4) 
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and 



q 2 {z)-qkz)^q-kz)=R{z)-(\-^) 
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or, alternatively, 



g(z)A 
G(z) 



G 2 W 



,R{z)-Q 2 {z) 



Q 2 (z) 



Introducing the complex variable 

C= [ Z Q(z)dz 



1 (6) 



(7) 



where, typically, {zq : Q 2 (zo) = 0}, one finds, after some cal- 
culations, that equation (6) can be rewritten 



Q(z) 



-x 2 



q(z)XY d 2 (q{z)iy 



Q 2 (z) 



Q{z) J dC 2 V Q(z) 
-I 2 — 



, 2 R(z)~Q 2 (z) j2 1 d 2 1 
A ^ , s A — 5 — — 



[cf. equation (2.2.5), with (2.2.1), in Ref. 4]. Using equa- 
tion (7), and the equality 



1 d 2 1 



1 d / 1 „ d 1 



q\ ( z ) dz 2 q\ ( z ) Q(z) dz\2 dz z s Q{z 



s(s-2) 

Y dzz s Q(z)J + 4z 2 2 2 (z) 



-iV 1 



(9) 



this can be rewritten 



Q(z) 



2 _ xZ (q{z)X\ i d 2 (q{z)kY* 



Q(z) J dZ 2 \Q(z 



■AM ^-^1-^-^^)1=1 do) 



G 2 « 



where 



, , 1 , d 1 
Ps(z) = =2* 



2 dz z s Q{z) 



(11) 



This is our new platform function. It allows us to write the 
right hand member of (3.5c) in Ref. 3 as 



R(z)-Q 2 (z) , dP ( Z/ d2 



Q 2 (z) dC 



-Po(z) 



(12) 



In practice, the choice of s in P s (z) more or less determines 
the base function Q(z) as we shall show next. 



III. BASE FUNCTION 

To obtain a solution of equation (10), we make the formal 
Ansatz 



q(z)X 
Q(z) 



In 



(13) 



TABLE I. Properties of the platform function P s (z), equation (11), 
for various choices of s. 



s Platform function Base function 



Typical z range 







1 d 1 

^Q(z) 
, l.d i 

-21 



1 d z 21 

2z 21 dz Q(z) 



Q 2 (z)=R(z)- 
Q 2 (z)=R(z)- 



4,2 

1(1+1) 



—°° < z < °° 
< Z < °° 
< z < °° 



and find that 

Y = l 



1 , ldP, G 2 (z)-*(z)- 

Yo = P H 

2 2 s ^ 2 dC 2g 2 (z) 



4z2 



(14a) 



(14b) 



Using equations (3) and (13) with X = 1, and truncating 
after « = 2, we obtain the approximate local solutions of the 
original differential equation ( 1 ) as 



V± = {Q(z)(l+Y 2 (z))}-?e^f Z Q^ l+Y ^ 



(15) 



It should be noted that the integrands of these solutions con- 
tain the total derivative 



2 G(z) ^T 



1 dP s (z) 

2 dz 



(16) 



which is integrable irrespective of the expression for R(z) in 
the original differential equation (1) and of the choice of the 
base function Q(z). Hence, we can write 



, ldP s / 1 2 \ Q 2 -R- 



2 dz 



s(s-2) 



2Q 2 



Q 2 -R- 



(17) 



From this we conclude that an obvious choice for the (square 
of the) base function is 



Q 2 (z)=R(z)- 



s(s-2) 
4z 2 



(18) 



For the new platform function (11), three cases can be dis- 
cerned (see also Table I): 

1. s = 0. This corresponds to the case of an unmodified 
base function, Q 2 (z) = R{z). One can show that the use 
of the platform function Pq(z) simplifies the calculation 
of phase-integral approximations, e.g., for the Weber 
functions. 

2. s = 1. This corresponds to the well-known Kramers- 
Langer modification for obtaining the base function 
g 2 (z) = R(z) - l/(4z 2 ). One can show that the use 
of the platform function Pi (z) simplifies the calculation 
of phase-integral approximations, e.g., for the Coulomb 
wave functions. 
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3. s = —21. This corresponds to the omission of the cen- 
trifugal barrier in the base function, i.e., the choice 
Q 2 ( z ) = R( z )+l(l + l)/z 2 . This is (for / ^ 0) a less 
known alternative; see, e.g., Chapter 7 in Ref. 3. One 
can show that the use of the platform function P_2/(z) 
simplifies the calculation of phase-integral approxima- 
tions for spherically and cylindrically symmetric prob- 
lems. 

In all the three above cases, the third-order solution becomes 

xe^+J^H^H (19) 

which exhibits the simplification attainable if the new plat- 
form function P s (z) is used. 

IV. CONCLUSIONS 

We have derived a platform function P s (z) for the third- 
order correction of the phase-integral method. Since P s (z) is 



readily generalized to problems involving more than one sin- 
gular point in R(z) and makes computational simplifications 
possible, this new platform function is in many cases prefer- 
able to the conventional PIM platform function Eo(z). 
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